)28 

1414 


'% 


9j. 


mi^-- 


■■**'f*»!»i(Wf 


WORKING  PAPER 
ALFRED  P.  SLOAN  SCHOOL  OF  MANAGEMENT 


A  Two  Priority  M/G/1  Queue  with  Feedback 

by  J.  Keilson 
GTE  Laboratories  Incorporated  and  Mass.  Inst,  of  Tech. 

and  L.D.  Servi 
GTE  Laboratories  Incorporated 

WP#1994-88 


MASSACHUSETTS 

INSTITUTE  OF  TECHNOLOGY 

50  MEMORIAL  DRIVE 

CAMBRIDGE,  MASSACHUSETTS  02139 


A  Two  Priority  M/G/1  Queue  with  Feedback 

by  J.  Keilson 
GTE  Laboratories  Incorporated  and  Mass.  Inst,  of  Tech. 

and  L.D.  Servi 
GTE  Laboratories  Incorporated 

WP#1994-88 


0^ 


\m- 


ABSTRACT 

An  M/G/1  queueing  system  with  tow  priority  classes  of  traffic  and  Bernoulli  feedback  for 
the  low  priority  class  is  studied.  For  the  low  priority  customers,  the  distributions  of  the 
ergodic  time  to  completion,  the  first  pass  time,  and  the  number  in  the  system  are  found  in 
the  transform  domain  and  the  first  moments  are  displayed.  The  distributed  Poisson 
approximation  introduced  previously  by  the  authors  is  them  employed  to  analyze  systems 
with  a  preempt-resume  high  priority  clocked  schedule  and  low  priority  traffic  with 
feedback. 
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1.  Introduction 

The  following  system  which  we  will  refer  to  as  the  feedback  system  is  studied:  a)  the  system 
consists  of  two  classes  of  traffic,  Cj  and  C2.  ^^ch  having  Poisson  arrivals  and  i.i.d.  service  times; 
b)  the  C]  customers  have  preempt-resume  priority  over  C2;  c)  upon  completion  of  service  of  any 
C2  task,  that  task  feeds  back  with  probability  9  to  the  end  of  the  C2  queue  to  again  await  service. 
This  system  contains  as  special  cases:  i)  the  M/G/1  system,  ii)  a  two  class  M/G/1  system  with 
prempt-resume  priority  as  considered,  for  example,  by  Keilson  and  Sumita  [9],  and  iii)  an  M/G/1 
system  with  feedback  considered,  for  example,  by  Takacs  [11]  and  Disney  [1],[2].  The  elements 
of  priority  and  feedback  interact  in  the  feedback  system  and  compound  the  effort  required  for  its 
analysis.  Many  of  the  formulae  obtained,  however,  can  be  given  a  structurally  simple  and 
informative  form. 

The  study  was  motivated  by  concern  for  the  performance  of  a  telecommunications  switch  with 
priority  and  feedback  featiires.  The  distributed  Poisson  approximation  introduced  in  Keilson  and 
Servi  [7]  and  the  results  obtained  for  our  feedback  system  permit  exact  analysis  of  a  system 
which  is  a  fair  approximation  to  one  such  telecommunications  switch  . 

In  Section  2,  notation  is  first  introduced.  Next,  using  analysis  of  a  system  without  feedback, 
the  p.g.f.  of  the  number  of  C2  customers  in  the  system  in  the  feedback  system  is  found  as  well  as 
the  mean  time  in  the  system  for  the  C2  customers.  Using  an  approach  related  to  that  of  Keilson  and 
Kooharian  [6]  the  equations  for  the  state- space  motion  are  found  in  Section  3  and  solved  in  Section 
4.  In  Section  5  this  analysis  is  used  to  find  the  transform  of  the  joint  distribution  of  the  duration  of 
the  first  pass  time  in  the  system  and  the  number  in  the  system  at  the  end  of  the  first  pass  time.  This 
leads  to  the  Laplace  transform  of  the  first  pass  time  and  its  expected  value.  In  Section  6  the 
transform  is  found  of  the  joint  distribution  of  the  sojourn  time  of  the  low  priority  customers  and  the 
number  in  the  system  at  the  epochs  at  which  low  priority  customers  leave  the  system.  In  Section  7 
the  busy  period  is  examined.  Finally,  in  Section  8  the  analysis  of  Section  2  is  applied  to  preempt- 
resume  clocked  schedules  with  feedback  of  low  priority  customers  using  the  distributed  Poisson 
approximation  introduced  analysis  in  [7].  The  results  are  then  validated  numerically. 

2.Notation.  The  mean  time  in  system  for  low  priority  customers. 

The  following  notation  will  be  helpful  for  reference.  Each  row  describes  an  M/G/1  system 
whose  Poisson  input  stream  is  the  class  listed  in  the  first  column. 
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service 

service  time 

ergodic 

Class 

arrival  rate 

time 

c.d.f 

LST 

utilization 

waiting  time 

Ci 

^l 

Ti 

Ayi(x) 

aj.j(s) 

Pi 

Wi 

C2 

h. 

T2 

At^(x) 

CXt^(s) 

P2 

W2 

C2EFF 

X2 

T2EFF 

At2EFf('^) 

047EFf(s) 

P2HFF-P2(1- 

-e)- 

■'      W2EFF 

Ci2 

X1  +  X2 

T12 

Ati2(x) 

0712(5) 

Pl2=Pl+P2 

W12 

Cs 

X1  +  X2 

Ts 

Ats(x) 

Ot-sCs) 

PS=Pl+P2(l- 

-e)- 

-1        Ws 

TABLE  2.1 

The  subscript  12  designates  the  priorityless  system  of  Cj  and  C2  customers  with  no 
feedback  .  The  corresponding  service  time  transform  is  given  by 

«T12(s)=7-V-«Tl(s)  +  .      \      «T2(s)  •  (2-1) 

1"''^2  "■l"^"'2 

The  subscript  2EFF  designates  a  class  without  feedback  whose  effective  service  time  is  a 
compound  geometric  mixture  of  multiple  convolutions  of  A^Cx)  with  itself.  The  corresponding 

service  time  transform  is  given  by 

aT^PPp(s)  =  (l-e)y  e*a'+V2(s)    =^^'^^"t2(s)  (2.2) 

T2EFFV  ^^  T2  l-Ba^^Cs) 

This  may  be  thought  of  as  arising  froin  the  class  of  customers  C2  feeding  back  to  the  head  of  the 
line  with  probability  0 . 

The  subscript  S  designates  the  priorityless  superposition  of  the  classes  Cj  and  C2EFF.  The 
corresponding  service  time  transforni  is  given  by 

«Ts(s)  =  r-V"Tl(s)  +  r-^aT2EFF(s)- 
A,i+A,2  "-I  '  ""2 

Equivalently,  from  (2.2), 

XiaTi(s)(l-eaT2(s))  +  X2(l-9)aT2(s)  ^. 

^^  (Xi  +  X2)(l-eaT2(s)) 
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This  should  be  contrasted  with  the  feedback  system  of  focal  interest  where  the  C2  customers 
return  to  the  back  of  the  line  with  probability  6  and  C^  has  preempt-resume  priority  over  C2 . 

The  follov/ing  representative  notation  will  also  be  employed. 

N(t)  the  number  of  C2  customers  in  the  system  at  time  t. 

T  the  ergodic  time  in  system  for  C2  customers  in  the  feedback 

system. 

9  the  probability  of  feedback  for  C2  customers. 

Tgpi  the  busy  period  for  Cj  customers  in  the  absence  of  Cj  ■ 

aBPi(s)  E[exp(-sT3pi)] ,  the  busy  period  transform  for  C^  customers. 

<t>PKi2(s)  E[exp(-sWj2)],   the  waiting  time  Pollaczek-Khinchin  transform 

for  the  priorityless  C12  system. 

4>PKs(s)  E[exp(-sWs)],   the  waiting  time  Pollaczek-Khinchin  transform 

for  the  priorityless  Cg  system. 

An  upper  case  asterisk  will  always  designate  forward  recurrence  time,  e.g. 
*  1  -  a-p,9(s)  *  1  -  a^-cCs) 

«T12(s)   =      E[Ti2]s  ^'^      "tS^s)   =      E[Ts]s  (2.4) 

With  this  notation 

<t>PKi2(s)= ^Ar and  (})pKs(s)  = '-A (2.5) 

1  -  pi2a^j2(s)  1  -  psa.p5(s) 

Case  1:  The  priority  system  with  no  feedback  (6  =  0) 

The  priority  system  with  no  feedback  has  been  treated  by  Keilson  and  Sumita  [9]  and  the 
distribution  of  the  low  priority  customer  delay  is  exhibited  there.  A  much  simpler  derivation  of 
this  distribution  is  given  next  which  sets  the  stage  for  the  more  difficult  feedback  case. 

The  arrival  of  a  tagged  C2  customer  in  steady  state  finds  an  ergodic  backlog  distribution  equal 
to  that  of  a  priorityless  system  of  Cj  and  C2  traffic  having  arrival  rates  Xj  and  X2  and  service 
times  c.d.f.  of  Ajj(x)  and  Ay2(x).  From  TABLE  2.1  the  L-S  transform  of  this  backlog  is 
<()p^j2(s)-  All  subsequent  C2  customers  have  no  effect  on  the  delay  time  of  the  tagged  C2 
customer.  However  subsequent  Cj  arrivals  modify  the  delay  time  as  interruptions  of  rate  Xi  with 
each  interruption  time  equal  to  that  of  Tbpi,  the  Cj  busy  period.  It  follows  at  once  that  the  C2 
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waiting  time  distribution  has  an  Laplace  transform  given  by  <t>pKi2(s  +  ^i  -^i<7bpi(s))  [4],  [5]. 
Using  similiar  reasoning  the  time  in  the  system  has  a  Laplace  Transform 

<h-i2(s)  =  <t>PKi2(s  +  ^1  -^i<7bpi(s))  ct^{s  +  ^j  -A^jObpiCs))  •  (2.6) 

Case  2:  The  feedback  system. 

The  process  N(t)  is  insensitive  to  a  certain  variant  on  the  underlying  order  of  service 

discipline.  This  insensitivity  enables  one  to  employ  Little's  Law  [10]  to  evaluate  the  expectation 
EfT^y^]  directiy. 

Theorem  2.1: 

a)  The  system  is  stable  when  p^  =  p  j  +  pj  (l-O)"""  <  1  and  the  idle  state  probability 
at  ergodicity,  £(»=),  is  given  by  1  -  p^. 

b)  The  probability  generating  function  of  N(oo)  is  given  by 

h(u)  =   (l)pKs(C(u))  aT2EFF(C(u))  (2.7) 

where 

C(u)  =  ?i2(l-u)  +  Xi[  1  -  aBpi(X2(l-u))  ],  (2.8) 

and  aj2EFF(^)  ^"^  ^KS^^^  ^^  defined  in  (2.2)  and  (2.5). 

c)  E[N(oo)]  =X2E[T5y3],  (2.9) 

and 

d)  For  a  low  priority  customer,  the  mean  time  in  the  system  is  given  by 

E[T     ]  =    (l-9)^lE[T,2]  +?i2E[T22]  +  2E[T2](1-Pi-P2) 

'^'  2(i-Pi){(i-e)(i-Pi)-P2} 

Proof:  The  ergodic  number  of  low  priority  customers  in  the  system,  N(«»),  would  be  the  same  in 
distribution  if  the  low  priority  customers  that  feed  back  were  to  feed  back  to  the  front  of  the  tine 
rather  than  to  the  back  of  the  line.  Hence,  h(u)  must  also  be  the  same  under  either  schedule. 
System  behavior  under  the  front  of  the  line  schedule  is  equivalent  to  that  for  a  two  priority  system 
in  which  the  low  priority  customers  have  an  arrival  rate  X,2  and  a  service  time  T2EFF  having  the 
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Laplace  transform  given  in  (2.2).  Using  the  reasoning  of  (2.6)  the  time  in  system  for  the  low 
priority  customers  of  a  two  priority  M/G/1  system  with  arrival  rates  X^  and  X2  and  service  times 
Tj  and  T2£pp  is  given  by 

E[exp(-sT5y5)]  =  <})pKs(s  ■•"  ^1  -^i^^BPi^^))  «t2EFf(s  ■•"  ^1  -^i<^Pi(s))  (2. 1 1) 

where  <l)pKs(s)  and  cjgpj(s)  are  defined  in  TABLE  2.L 

From  the  distributional  form  of  Little's  Law  [8],  provided  for  conveniencein  Appendix  A.l,  the 
distribution  of  N(o°)  and  T^     are  related  to  each  other  according  to 

h(u)  =E[uN(-)]  =  E[exp(-(X2(l-u))T3y,)]  (2. 12) 

so  (2.7)  follows. 

Equation  (2.9)  follows  from  either  Little's  Law  [10]  or  from  differentiation  of  (2. 12).  From  (2. 1 1) 

H[T3,]    =B^^^dl^^l2^  (2.13) 

1  -pl 
where  W5  ,  T2Epp  and  pi  are  defined  in  TABLE  2.1.  Equation  (2.10)  then  follows  from  (2.13) 

and  the  Pollaczek  -  Khintchine  formula, 

E[WS1  =  -^  1^  (2.14) 

where  Ts  and  ps  are  defined  in  TABLE  2.1.  (Alternatively  (2.10)  follows  from  (2.7)  and  the 
observation  that  E[N(oo)]  =  h'(l).)* 


Remark:  Note  that  when  Xj  =  0  and  9=0,  (2.10)  coincides  with  E[Tsys]  for  the  ordinary  M/G/1 
queue.  When  X^  =  0,  (2. 10)  agrees  with  E[Tsys  ]  for  the  M/G/1  queue  widi  feedback  as  given  by 
[11,  equation  (35)].   When   6=0,  agreement  is  found  with  E[Tsys  ]    for  the  preempt-resume 

M/G/1  queue  with  two  classes  of  traffic  given  in  [9,  equation  (3.8f)]. 

In  the  case  of  no  high  priority  customers,  i.e.,  Xj  =  0,  equation  (2.7)  is  equivalent  to 

(l-u)(l-e-p2)  0x2^2(1 -u)) 


h(u)  = 


euaj2(^20-"))  +  (l-6)aT2(^2(l-"))  -  " 


so 


-(l-P2(l-e)-0u(l-aT2(M^-"))) 
h(u)  -  E(c«)    = 
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which  is  consistent  with  [11,  equation  (20)] . 

3-    The  motion  of  the  system  on  its  state  space 

To  describe  the  feedback  system,  the  following  notation  will  be  employed. 

Let  N(t)  be  the  number  of  low  priority  customers  in  the  system. 

Let  X(t)  be  the  virtual  time  to  availability  of  service  for  C2  customers  who  are  not  in  service 
at  time  t,  i.e.  let 

X(t)  =  Bi(t)  +  R2(t), 
where  B  j(t)  is  the  C^  (high  priority)  backlog  at  time  t  and  R2(t)  is  the  residual  service  time  of  any 
C2  customer  who  has  been  interrupted  and  is  awaiting  resumption  of  service. 

The  reader  will  note  that  [N(t),  X(t)]  is  not  bivarate  markov  for  the  following  reason.  When 
X(t)  becomes  zero  at  a  C2  service  completion  epoch,  N(t)  decreases  by  one  .  On  the  other  hand, 
when  X(t)  becomes  zero  at  a  Cj  service  completion  epoch  (when  no  C2  customers  having  been 
served  in  the  current  busy  period),  N(t)  does  not  change.  A  third  process  J(t)  taking  two  values  S 
and  V  ,  respectively,  distinguishes  between  such  system  histories.  Specifically: 

a)  J(t)  =  V  if  the  server  is  busy  and,  in  the  current  busy  period,  no  C2  service  has  been 
provided; 

b)  J(t)  =  S  if  the  server  is  busy  and  some  C2  service  has  been  provided  during  the  current  busy 
period. 

Finally,  let  E  be  the  idle  state. 

With  this  notation  it  is  clear  that  [N(t),  X(t),  J(t)]  is  a  multivariate  Markov  process  on  the  state 
space  N  =  V  +  S  +E    with 
S  ={(n,x,S) ;  1  <  n  <oo,  0  <  X  <  «};  V  =  {(n,x,V);0  <  n  <oo,  0  <  x  <  00);  E  ={E}. 

Let  Tj  be  the  service  time  of  a  type  i  customer  having  density  ajj(x),  i  =  1,2  Jumps  on  the 
state  space  N  induced  by  arrivals,  busy  period  terminations  and  low  priority  service  completions  are 
enumerated  in  TABLE  3.1. 
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Epoch  Changes    from  V  or  E  Changes  from  S 

Cj  arrival  (n,x,V)  -->    (n,  x  +  Tj,  V);  (n,x,S)  -->  (n,  x  +  Tj,  S); 

E  ->    (l.Ti.V) 

Cj  arrival  (n,x,V)  -->  (n+1,  x,  V)  ;  (n,x,S)  -->  (n+1,  x,  S); 

E     ->(1,T2,S) 

termination  of  the        (n,0,V) -->  (n,T2,S) 
J(t)  =  V  mode 

C2  service  completion  (n,0,S)  -->  (n-l,T2,S) ;  n  >  1 

with  no  feedback  (1 ,0,S)  -->  E 

C2  service  completion  (n,0,S)  -->  (n,T2,S) 

v.'iih  feedback 

TABLE  3.1 

When  underlying  service  time  distributions  are  absolutely  continuous  the  distribution  on  the  state 
space  A^  can  be  described  in  terms  of  densities  on  S  and  V.    Let 

fvn(x,t)  =  ^Prob[  N(t)  =  n,  X(t)  <  x  ,  J(t)  =  V] ; 


E(t)  =  Prob  [  System  is  in  the  idle  state  E  at  time  t] . 

The  equations  of  motion  on  the  state-space  may  be  written  down  in  the  customary  way  (see,  e.g., 
[6]). 

One  finds  that 

^  =  -  (?Li+?L2)  E(t)  +  fvo(0,t)  +  (l-e)fsi(0,t)  (3.1) 
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^vo(x,t)      3fvo(x,t) 
3i  5x 


=  -^.jfvo^'^'O  +  ^1^(0  aTi(x)  -  ^ifvo(x.t)  +  Xifvo(x,t)*aTi(x),  (3.2) 


^fvn(x.t)       3fvn(x.t) 

5t  5x 


= -X2fvn(x,t)  +  >^fvn.i(x,t)    -  Xifv„(x.t)  +  Xifvn(x,t)*aT-i(x),     n^l,  (3.3) 


afg  (x,t)    afs  (x.t) 

5t  5x 


=  -  Vsn^'^'t)  +  Vsn-l(X'0    -  ^lfsn(X'0  +  >.ifsn(x,t)*aT.i(x) 

+  fvn(0,t)  aT^(x)  +  (l-0)fs„^i(O,t)aT^(x)  +  efs„(0,t)aT2(x),    n  >  1  (3.4) 

^fsi(x.t)      3fsi(x,t) 
5t  oix 

=  -X2fsi(x,t)  +  X2E(t)aTV2(x)    -Xifsi(x,t)  +  Xifsi(x,t)*aji(x) 

+  fvi(0,t)aT2(x)  +  (l-e)fs2(0,t)a^(x)  +  efsi(0,t)aT^(x)  (3.5) 

and 

E(0)  =  1.  (3.6) 

Equation  (3.1)  for  example  states  that  the  rate  of  change  in  E(t)  is  due  to  losses  associated  with 
the  Cj  and  C2  arrivals  and  gains  from  busy  period  terminations.  Equation  (3.4)  says  that  the  time 
rate  of  change  of  fY„(x,t)  seen  by  an  observer  moving  with  velocity  - 1  consists  of  losses  from 
{(n,x,S)}  due  to  arrivals  of  C2  customers,  gains  from  {(n-l,x,S}  due  to  C2  arrivals,  internal 
losses  and  gains  on{(n,x,S)  due  to  Cj  arrivals,  gains  from  (n,0,V)  for  termination  of  busy 
periods  initiated  by  a  Cj  customer  arrival,  and  service  completion  of  C2  customers.  The  other 
equations  have  a  similar  probabilistic  interpretation. 

4.  Ergodic  analysis 

In  the  next  theorem  the  transform  of  the  joint  distribution  of  N(o«),  the  number  of  low  priority 
customers  in  the  system  at  ergodicity,  and  of  X(oo)  are  obtained  from  Equations  (3.1)-(3.6). 
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Let 

oo  oo 

(I>Voo(u,x)  =  X  fvn(x,-)  u"  ;    <I>Voo(u.w)   =  L  [  5;  fvn('^.~)  ""] 
n=0  n=0 

00  oe 

(t)s^(u,x)  =  £  fs„(x,oo)  u"  ;    Os^(u.w)  =  L  [  ^  fj^Cx.-)  u"] 
n=0  n=0 

(where  L[]  is  the  Laplace  Transform  operator)  and  let   ^(u)  be  the  solution  of  the  functional 

equation 

C(u)  =  ?i2[l-u]  +  >.i  [1  -  OtiCCCu))  ] .  (4.1) 

This  functional  equation  plays  a  key  role  in  the  solution.  As  shown  in  Appendix  A.2,  the  functional 
equation  (4.1)  has  the  unique  solution  given  in  (2.8). 

Theorem  4.1  (The  Ergodic  Distribution  on  the  State  Space) 

The  ergodic  distribution  of  [N(t),  X(t),  J(t)]  on  its  state  space  is  described  by 

^     ..  ^  u  E(oo)  C(u)  (  aT2(w)  -  aT2(C(u))  ) 

^~    '  [w->.2(l-u)-Xi(l-aTi(w))][u-a-r2(C(u))  +  eaT2(C(u))(l-u)] 


^     I      ^     1  c.  aji(C(u))-aTi(w) 

<I>V~(u.w)  =  ?.iE(oo) ■    ]\     ^     ■    ^/' -— -   .  (4.3) 

and 

E(oo)=l-ps=l-pi-p2(l-e)-l.  (4.4) 

Proof:  Equations  (4.2)-(4.4)  are  obtained  from  the  equations  of  motion,  (3.1)  -  (3.6),  by  standard 
methods  employing  p.g.f.'s,  Laplace  transforms,  regularity  conditions  and  algebra.  Details  are 
provided  in  Appendix  A.3.  ♦ 

Remark:  By  definition  h(u)  =  <I>5„(u,0)  +  <I>y^(u,0)  +  E(oo).  Hence,  from  equations  (4.2)-(4.4) 
and  some  algebra,  equation  (2.7)  follows.  Details  are  provided  in  Appendix  A.4. 
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5.  The  first  pass  time 

We  define  the  first  pass  time  of  a  Cj  customer  to  be  the  time  from  the  arrival  of  that  customer  to  the 
system  until  the  completion  of  that  customer's  first  service.  The  expected  value  of  the  first  pass  time 
at  ergodicity  is  a  measure  of  the  performance  of  the  system.  Its  distribution  is  also  an  important 
ingredient  for  the  calculation  of  the  total  time  in  the  system- 
Theorem  5.1  (Ergodic  First  Pass  Time  and  the  Number  in  the  System):  The 
transform  of  the  joint  distribution  of  V. ,  the  duration  of  the  first  pass  time,  and  of  N,,  the  number 

of  C2  customers  in  the  system  at  the  end  of  the  first  pass  is  given  by 
T^i(u,s)  =  E[u^l*  e-sVl] 

=  ot^^Cw)  [  <I>voo(z.w)  +  E(oo)]  +  <Ds^(z,w) .  (5.1) 


z=  z(u,s)  =  (l-e+Gu)  aiv2(w(u,s))  (5.2) 

w  =  w(u,s)  =  s  +  ^2(l-u)  +  X\-Xi  CTbpj(s  +  X.2  (1-u) ).  (5.3) 


Here 

and 

Hence, 

i&l(u,s) 

(l-(pj+p2))aj2(w)               irc/   N       6p2  . ,  X    (u-z)    ,,,,  ,  ,,     ,,  -, 

=  [  *^ — ]  [E(~)  H L_h(z)-4p^]/(l-(p  +P2))     (5.4) 

l-(Pi+P2)a.j,j2(^)+e(l-u)aT^(w)/w  (1-e)         ^^l^2J 

* 
where  a    ^(w)  and  h(z)  are  defined  in  (2.4)  and  (2.6). 

Proof:  Equation  (5.1)  follows  from  a  generalization  of  the  argument  of  Theorem  2.1.  Equation 
(5.4)  follows  from  algebra  involving  (4.2),  (4.3),  and  (5.1).  Details  are  provided  in  Appendix  A.5. 

The  Laplace  transform  of  the  distribution  of  the  first  pass  time  at  ergodicity  can  be  found  by  simply 
evaluating  (5.4)  at  u=l.  This  transform  is  given  in  the  next  corollary: 

Corollary  5.2. (Ergodic  First  Pass  Time  ) 

The  Laplace  transform  of  the  ergodic  first  pass  time  distribution  is  given  by 


page  11 


=    <l>pKi2(^(s))  aT^(w(s))  [p  +  (1-p)  h(a^(w(s)))  a^(w(s))  ]  (5.5) 

where 

w(s)  =w(l,s)=  s+X,i-XiaBpi(s), 

p  =  E(oo)  /  (I-P1-P2)  =  (l-pi-p2(l-e)-i)  /  (I-P1-P2), 

and  <}>pKi2(^)  ^^  given  in  (2.5). 

The  mean  is  given  by 

E[Wi]  =  -4^'i(0) 

^((l-e)(l-pi)-ep2){  X2E[T22]+(l-e)?ciE[Ti2]-2p2E[T2]  )    ^       E[T2](l-9)        ^^  ^^ 

2(i-pi)2(i-e)  ((i-e)(i-pi)-p2)  (i-e)(i-pi)-p2 

Proof;  Equations  (5.5)  and  (5.6)  follow  from  (5.4).* 
Corollary  5.3  (Joint  Distribution  under  no  feedback) 

ife=o 

i^i(u,s)  =  <})pKi2(w(u,s))    aj2(w(u,s))  (5.7) 

where  w(u,s)  and  <t)pKi2(^)  ^^  defined  in  (5.3)  and  (5.6). 

Proof:  Equation  (5.7)  follows  from  (5.4).* 

Remark:  In  thfe  case  of  no  feedback  ,  i.e.,  6=0,  from  (5.7)  and  (2.6) 

»Fi(s)=  i3i(l,s)  =  <lvri2(s)  (5.8) 

which  is  equivalent  to  equation  (3.5)  and  (3.7)  of  [9].  K  6  =  0  and  X\=  0  this  equation  simplifies  to 
*Fj(s)  =   <)>pKi(^)  01^2(5)  as  expected. 

Littie's  Law,  [10],  can  be  used  to  relate  the  average  number  of  customers  in  a  system  to  the  average 
time  in  the  system.  The  following  corollary  verifies  the  distributed  form  of  this  law  for  the  case  of 
9=0. 
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Corollary  5.3  (A  Distributional  form  of  Little's  Law): 
K  e  =  0.  then 


Hence 
and 


'¥y(k2  -  >.2U)  =  Oi(l,>.2  -  ^2u)  =  t5i(u,0)  =  h(u).  (5.9) 

E[N(oo)]  =  >.2E[Tsys]  (5.10) 

Var[T^3]  =  Var[X2N(oo)]  +  E[N(oo)].  (5.11) 


Proof:  If  e  =  0  then,  from  (2.8),  (5.3)  and  (5.7), 

^l(X2-X2u)=  l5i(l,X2->^2u)=l»i(u,0)=    <t)pKi2(C(u))   O^iW)- 
But  if  0  =  0  then  cXy2EFF(s)  =  aj2(s),  <t>pKi2(^)~  ^KS^^^  ^°  ^°™  (^-^^  ^^  equals  h(u)  so  (5.9) 
follows.  If  0  =  0  there  is  no  feedback  so  the  time  in  the  system,  Tsys,  is  equal  to  the  first  pass  time, 
V*.  Therefore  "¥{(0)  =  E[Tsys],  ^i"(0)  =  E[Tsys2],  hXD  =  E[N(oo)],  and   h"(l)  =  E[N(~)2]  - 

E[N(oo)].  Hence  (5.10)  and  (5.11)  follow  from  (5.9).* 
6.  Total  time  in  system 

Theorem  6.1:  Let  x(u,s)  be  the  transform  of  the  joint  distribution  of:  a)  the  number  of  customers 
remaining  in  the  system  at  departure  and  b)  the  total  elapsed  time  from  arrival  to  the  system  to 
departure  from  the  system.  Then  x(u,s)  satisfies 

X(u,s)  =  (l-0)i^i(u,s)  +  0aT-2(w(u,s))  x(z(u,s),w(u,s)))  (6.1) 

where  z(u,s),  w(u,s)  and  t3i(u,s)  are  defined  in  equations  (5. 2), (5. 3)  and  (5.4) . 

Proof:  Using  an  argument  similar  to  Theorem  (5.1)  one  can  show  that  'dj^.i(u,s)  = 
otp2(w(u,s))i5j(z(u,s),w(u,s))  where  i^j(u,s)  is  the  joint  transform  of  the  duration  of  the  first  i 

passes  and  the  number  of  C2  customers  in  the  system  at  the  end  of  the  ith  pass.  But    x(U'S)  = 
00 

^  (l-0)0'"^dj(u,s)    so  equation  (6.1)  follows.  Details  are  provided  in  Appendix  A.6.* 
i=l 


7.  The  busy  period 
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The  following  theorem  gives  a  generalization  of  the  Takacs  busy  period  equation. 
Theorem  7.1:  (The  Busy  period)The  busy  period  has  a  Laplace  transform 

<^BP(S)  =  : ^Cy*BPl(s)  +   ,  .      <7*BP2(S)  (7.1) 

1  ''"      2  1  "''      2 

a*BPi(s)  =  aT.i[  s  +  X,  -  >.jCr*Bpi(s)+  X^  -  X^^s^^^^is)  ],  (7.2) 

<^BP2(S)  =  0^p2EFf[  S  +  ^1  -  XjO*gpi(s)+  X^  '  ^2'^BP2(s)  1'  C^-^) 


where 


a^j(s)  and  aj2EFF(s)  are  defined  in  TABLE  2.1, 
and 

E[T3,i=  p.^P2('-e)-'     _1_^EI±  ,,.4, 

l-Pi-P2(l-9)-i     X,U2       1-Ps 
Proof:  See  Appendix  A.7. 

Note  that  if  Xj  or  Xj  is  set  equal  to  zero  then  equations  (7.2)  and  (7.3)  reduce  to  the  familiar  Takacs 
busy  period  equation. 


8.   Clocked  schedules 

Consider  a  system  in  which  the  priority  tasks  arrive  with  deterministic  interarrival  times  kA. 
Suppose  further  that:  a)  the  service  time  of  the  priority  tasks,  Ti  is  infinitely  divisible,  i.e.,  for 
every  integer  N,  Ti  is  distributed  as  the  sum  of  N  independent  and  identically  distributed  random 
variables;  b)  the  low  priority  tasks  arrive  with  an  exponentially  distributed  interarrival  time  and 
have  a  service  time  with  a  general  distribution;  c)  the  low  priority  tasks  feed  back  with  probability 

e. 

The  distributed  Poisson  approximation  of  [7]  may  be  used  to  model  the  high  priority  traffic  as  a 
sequence  of  Poisson  streams  parametrized  by  the  variable  N  with  rate  ^in  =  N/A  and  service  time 
having  Laplace  Transform  a  in  (s)  =  aTi^'^(s). 
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As  in  the  classical  proof  of  de  Finetti's  Theorem  [3],  one  has 

aTi(s)  =  limj,^_^„[exp  {  -  >.in  A  [1  -  am  (s)  }  ],  (8.1) 

i.e.  the  distribution  of  the  high  priority  tasks  arriving  in  A  time  units  under  the  approximate 
Poisson  model  with  arrival  rates  ^in  and  service  times  aiisKs)  is  exactly  ajiCs).  From  (8.1)  it  can 
be  shown  that 

limN_^^{>.lNAE[TiN]}  =  E[Ti],  (8.2) 

lim  N_,^  {  XiN  A  E[Tin2]  }  =  Var[Ti],  (8.3) 


and 


Pl  =  lim  N^oo  {  ^'IN  E[Tin]  }  =  E[Ti]  /  A.  (8.4) 


Equation  (8.1)  suggests  that  the  clocked  schedule  with  feedback  could  be  accurately  modeled  as  the 
limit,  in  N,  of  a  system  with  two  classes  of  Poisson  traffic  having  arrival  rates  ^in  and  X2,  service 
times  with  Laplace  transforms  aiN(s)  and  a2(s)  and  having  a  feedback  parameter  9.  In  fact,  in  [7] 
this  distributed  Poisson  approximation  has  been  shown  analytically  as  well  as  numerically  to  be 
highly  accurate  for  the  case  of  6=0.  The  mean  time  to  completion  for  the  limiting  system  can  now 
be  readily  obtained: 

Theorem  8.1  (Mean  Time  to  Completion  of  the  Distributed  Poisson  Approximation 
to  the  Clocked  Schedule  with  Feedback) 

The  mean  time  to  completion  is  given  by 

(l-e)Var[Ti]/A  ^X^EIT^^]  +  2E[T2](1-Pi-P2)  ,.  ^, 

'^'  2(i-Pi){(i-e)(i-Pi)-P2) 

Proof:  Equation  (8.5)  follows  immediately  from  equations  (2.10),  (8.2),  (8.3)  and  (8.4).  ♦ 
Results 
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To  test  the  accuracy  of  (8.5)  the  following  three  examples  were  compared  with  simulations: 

CASE  1:  X-i  =  1;  Ti  has  an  Erlang-3  distribution  with  a  mean  of  .2;  Xj  =  1;  T2  has  an  exponential 
distribution  with  a  mean  of  .15;  Four  values  of  6  were  examined:  9=0,. 25, .5,  and  .75. 

CASE  2:  X,i  =  1;  Ti  has  an  Eriang-2  distribution  with  a  mean  of  .2;  X2  =  5;  T2  has  an  Erlang-2 
distribution  with  a  mean  of  .05;  Five  values  of  9  were  examined:  0=0,  .1,  .2,  .3,  .4,  and  .5. 

CASE  3:  Xi  =  1;  Ti  has  an  exponential  distribution  with  a  mean  of  .2;  X2  =  .1;  T2  has  an  Eriang-5 
distribution  with  a  mean  of  .2;  Eight  values  of  9  were  examined:  9=0,  .25,  .50,  .75,  .90,  .95,  .96, 
and  .97. 

In  TABLE  8.1,  for  each  case  the  value  of  E[T     ]  is  computed  via  simulations.  The  simulated  value 

of  EnTgys],  E[Tg  g]  .  ,  is  given  with  a  95%  confidence  interval.  In  addition,  the  theoretical  value 

of  E|T     ]  from  (8.5),  E[T     ]^^   ,  is  given  along  with  the  relative  error, 
IE[T,^J  .   -  E[T,^J,.       I       ^"^ 

*-    sys-'sim       '^    sys^ theory 


EtTsysJsin. 
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CASE         e  E[T,y,l,,^^^y 


^^^sysJslm  ETT— T~^ 

•^     sys^slm 


1 
1 
1 
1 

2 
2 
2 
2 
2 
2 

3 
3 
3 
3 
3 
3 
3 
3 


00 

.244 

,25 

.347 

50 

.617 

75 

3.042 

00 

.107 

10 

.122 

20 

.144 

,30 

.177 

,40 

.233 

,50 

.349 

,00 

.284 

,25 

.374 

,50 

.554 

,75 

1.135 

,90 

3.342 

,95 

9.962 

,96 

16.583 

,97 

49.688 

.244±  .090 

.352±  .379 

.620±  1.060 

3.066±  7.251 

.108+  .251 
.122±.212 
.154±  .095 
.187±  .000 
.243±  .095 
.357±  .402 

.283±       .232 

.368±       .474 

.55 1±       .943 

1.1 28±     2.496 

3.327±     8.425 

9.962±    26.320 

16.232±  43.390 

50.171±136.89 

TABLE  8.1 


0.0% 
1.4% 
0.5% 
0.8% 

0.9% 
0.0% 

6.5% 
5.3% 
4.2% 
2.2% 

0.4% 
1.6% 
0.5% 
0.6% 
0.5% 
0.0% 
2.2% 
1.0% 


APPENDIX 


A.  1:  Theorem  (Keilson  and  Servi  [8]): 


Let  an  ergodic  queueing  system  be  such  that  for  a  given  class  C  of  customers, 


a)  arrivals  are  Poisson  of  rate  X2 


page  17 


b)  all  arriving  customers  enter  the  system,  and  remain  in  the  system  until  served 

c)  the  customers  leave  the  system  one  at  a  time  in  order  of  arrival 

d)  arriving  customers  do  not  affect  the  time  in  the  system  of  previous  customers 

then  N(oo)  and  Tgys,  the  ergodic  number  in  the  system  and  the  time  in  the  system,  are  related 
according  to  equation  (2.12). 

A.2  Proof  that  equation  (2.8)  is  a  solution  to  equation  (4.1): 

Let  Q(\)  =  V  +  Xi(l  -  C7gp2(v))  and  V(w)  =  w  -  ^i(l  -  ajj(w))  where  cy3pj(v)  is  the  Laplace 
transform  of  a  busy  period  of  an  M/G/1  queue  with  arrival  rate  Xi  and  service  time  Laplace 
transform  ajj(w).    It  can  be  easily  shown  that  Q(v)  and  V(w)  are  inverses  on  the  appropriate 

domain,  i.e, 

V(Q(v))  =  ri(v)-?ii(l-aT^i(Q(v))} 

=  v  +  ?ii  {  l-OBpi(v)}->.i{l-aTi(f2(v))}   =  V 
because  the  Takacs  equation  is  equivalent  to  cJbp2(v)  =  ajj(Q(v)). 

Hence,  on  the  appropriate  domain 
w  =  Q(V(w)) 
=  V(w)  +  ?ii(l-aBpi(V(w))) 
=  w  -  Xi  ( 1  -  a^i(w))  +  Xj  (  1  -  aBpi(w  -  Xi  ( 1  -  aT-i(w)) ). 

This  implies  that 

aji(w)    =  agpi(w  -  Xi  (  1  -  aTi(w)) ).  (A.l) 

Evaluating  (A.l)  at  w  =  ^(u)  gives 

aTi(C(u))   =OBpi(C(u)-)ii[l-aTi(C(u))]  ) 
so,  from  (4.1), 

aTi(C(u))  =  aBpi(?t2[l-u]).  (A.2) 

Therefore  equation  (2.8)  follows  from  (4.1)  and  (A.2).  ♦ 

A.3  Proof  of  Theorem  4.1: 

From  equations  (3.3)  and  (3.4)  the  Laplace  transform  with  respect  to  the  variable  x  and  the 
probability  generating  function  with  respect  to  the  variable  n  is  found  and  its  limit  is  evaluated  when 
t  approaches  <«.  Then 

<)v^(u,  0)  -  w  Oy<^(U,W  ) 

=  -X2(l-u)Ovoo(u.w)  -XiOy^(u,w)  +  XiOv^(u,w)aTi(w)  +  XiE(co)aji(w), 
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so 


({V^(u,0)  -  >.iE(oo)aTi(w) 

^eo(u,w)  = -—  (A.3) 

w  -  /.2(1"U)  -  ^i(l  -  aji(w)) 

From  (4.1)  the  denominator  of  (A.3)  vanishes  at  w  =  ^(u).  Regularity  in  w  in  the  half-plane  {w: 
Re(w)  >0  )  requires  that  the  numerator  also  vanish  there.  One  then  has  from  equation  (A.3) 

(>V^(u,  0)  =  >.iE(oo)  aTi(C(u)).  (A.4) 

Equation  (4.3)  then  follows  from  (A.3)  and  (A.4). 

dE(oo) 
Also,  from  (3.1),       }       =  0  .  Hence 

(Xl  +  >.2)E(~)  =    fvo(0.~)  +  (l-e)fsi(0,oo) .  (A.5) 

Equations  (3.4)  and  (3.5)  imply  that 

-  3(t)c^(u,x) 
=  A2(l-u)Os  Ju,w)  -  Xi(l  -  aTi(w))  Og  Ju.w) 

+  [  <tVoo(u.O)  -  fvo(0'~)  i  aT2(w)  +  ^2E(~)aT2(w)u 

+  (l-0)[  <l)soo(u,0)u-i  -  fs,(0,co)  ]aj2M  +  e(l)s^(u,0)aT2(w).  (A.6) 

From  (A.4),  (A.5)  and  (A.6)  one  obtains, 

^     ,      ,       <t>s«(u.O)  [  l-eaT2(w)-(l-9)u-iaT2(w)  ] 

<Pr^(u,W)   =      — ^ — 

^  w-X2(l-u)-Xi(l-aTj(w)) 


E(oo)aT2(w)  [  ^,(1-  Obpi(>>2(1-")))  +>-2(^-")  ] 
w-X2(l-u) -Xi(l-aTi(w)) 
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(A.7) 


From  (4.1)  the  denominator  of  (A. 7)  vanishes  at  w  =  ^(u).  Regularity  in  w  in  the  half-plane 
{w:Re(w)  >0  )  requires  that  the  numerator  also  vanish  there.  One  then  has  from  equation  (A.7) 

0  =  <t)soo(u,0)  [  1  -  eaT^(C(u))  -  (l-e)u-laT2(C(u))  ] 

+  E(oo)a^(^(u))  [  X^(  l-aBpi(X2(l-u)))  +X2(l-u)  ] 

which,  from  (2.7),  implies 

^     r    m  -  E(oo)aT.2(C(u))  C(u) 

(bo^(u,0)    = — : :  (A. 8) 

"^^^  1  -  eaT2(C(u))  -  (l-e)u-iaT2(C(u)) 


Equation  (4.2)  then  follows  from  (2.8),  (A.7)  and  (A. 8).  Equation  (4.4)  follows  from  equations 
(4.2)  and  (4.3)  and  the  condition  that  <I)s<^(l,0)  +  OycoCLO)  +  £(«>)  =  1.  ♦ 

A.4  Derivation  of  h(u)  from  equation  (4.2)-(4.4): 
From  (4.2)  and  (4.3) 

O3^(u,0)  =  ^ -E(oo)  C(u)  (1  -  aT2(C(u))  ) 


X2(l-u)[  1  -  eaT2(C(u))  -  (l-e)u-iar2(C(u))  1 
and 

^       .    m       T..    N           ^r    ,       Xi(l-aTi(C(u)))  +  Ml-") 
<I>V<«(u>0)  +  E(oo)    =     E(oo) i ti— ^ ^ 

/^(l-u) 
Then  from  (4.1) 

E(~)C(u) 


<Dvoo(u,0)  +  E(co)  = 


But 

h(u)  =  Os^(u,0)  +  Ov^(u,0)  +  E(oo) 

so,  after  some  simplification, 

h(u)    -  -E(oo)C(u)(l-9)aT2(C(u))  (^  9^ 

X2(  u  -  euax2(C(u))  -  (l-e)aT2(C(u))  ) 

-E(oo)C(u)(l-9)aT2(C(u)) 


X2(l-aT2(C(u)))  ->^(l-u)  (  l-eaT.2(C(u))  ) 
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which,  from  (4.1), 


which,  from  (2.3), 


-E(oo)C(u)(l-e)aT2(C(u)) 


>.2(l-aT2(C(u)))  +[?ti(l-aTi(C(u)))  -  C(u)  ]  (  l-eaT2(C(u)) ) 

E(oo)C(u) (l-e)aT^(C(u)) 

C(u)  -  (Xi+>^)(1-  (Xts(C(u)))        l-eaT^(C(u)) 


^  l^Ps (l-e)a-i^(C(u)) 

^  1  -  aTs(C(u))      l-8aT^(C(u)) 

"  ^^        C(u)E[Ts] 

equals  (2.7).  ♦ 

A.5  Proof  of  Theorem  5.1: 

The  joint  transform  of  V*,  the  duration  of  the  first  pass,  and  N*,  the  number  of  C2  customers  in 

the  system  at  the  end  of  the  first  wait  conditioned  on  the  value  of  (N,X,J)  upon  arrival,  if  no  future 
Cj  or  C2  customers  arrive  ,  is 

E[u^Vs^*  I  N=n,  X  =  X,  J  =  V]    =  (1-e+eu)"  aT2"+Hs)  e-"  .  (A.IO) 

The  first  term  reflects  the  fact  the  each  of  the  n  C2  customers  will  feed  back  with  probability  9. 

If  J  =  S  the  right  hand  side  of  (A.  10)  is  modified  to  reflect  the  fact  that  a  C2  customer  is  in  service 

so 

ELu'^VsV*  I  N=n,  X  =  X,  J  =  S]    =(l-e+eu)"aT2"(s)  e-",  (A.ll) 

and 

E[u^Vs^*  I  system  is  idle  i.e.,  N(oo)=0,  X  =  0,  J  =  E]    =  a-n(s).  (A.12) 

Hence, 


00 
00 


E[u^VsV*]=   X        f    (l-0+0u)"aT2""'ks)e-"  fvn(x,~)dx 
n=0   s=0 


00 
00 


+  I      f    (1-0+eu)"  aT2"(s)  e-"  fsn(x,«>)dx   +  aT2(s)E(-o).  (A.13) 

n=Os=0 
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Therefore, 

E[u^VsV*]  =  a^(s)  (Dv^[(l-e+eu)aT2(s),s]  -KDs^[(l-e+eu)aT^(s),s]  +  aT^(s)E(oo)  (A.  14) 


In  the  presence  of  future  Cj  arrivals,  the  number  of  C2  customers  in  the  system  at  the  end  of 
the  first  wait,  Np,,  is  still  N*.  However,  the  Ci  arrivals  modify  V   by  interruptions  of  rate  Ki  with 

each  interruption  time  is  equal  in  distribution  to  that  of  Tbpi,  the  Ci  busy  period.  Hence,  in  the 

presence  of  future  Ci  arrivals  the  joint  wait  transforms  of,  V, ,  the  time  to  completion  of  the  C2 

* 
customers  and  Npj  is  [4], [5] 

E[u^Pi*e-s^i*]  =  E[  u^*  e<^  +  ^1  -^i«^bpi(s))V*  ]  (A.15) 

where  CJgpj(s)  is  defined  in  TABLE  2.1. 


In  the  presence  of  both  future  Cj  arrivals  and  future  C2  arrivals  the  time  to  completion  C2 
customers,  Vi,  is  still  V..  However,  the  number  of  C2  customers  in  the  system  at  the  end  of  the 

first  wait,  Npi,  is  a  modification  of  Np.  which  reflects  the  Poisson  stream  with  rate  ^2  that  arrives 

* 
during  V, .  Hence  the  joint  distribution  of  Npi  and  Vj  is 

i^l(u,s)  =    E[u^Pi  e-sVi]  =  E[u^Pi*e-  (^  +  ^  ->^u)Vi*  j  (^  16) 

Combining  (A.  14),  (A.15)  and  (A.  16)  gives  (5.1)-(5.3). 

Fix)m  (2.8)  and  (4.3), 

.    ^r^       .        X      T..    M  aT2(w)E(°o)(w  -  C(z)) 

o^(w)  [ Oy^(z,w)  +  E(~)]  =  '"■  — - .  (A.  17) 

w  -  ^2(1-2)  -  Xi(l  -  aTi(w)) 

From  (4.2)  and  (A.  17), 
<^Soo(z.w)  +  aT2(w)  [  Oy^(z,w)  +  E(oo)] 

=  [ IM , 

w-X2(l-z)-Xj(l-a-p2(w)) 
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^      z  C(z)  (aT2(w)  -  aT2(C(z)))  ,   ,,       ^,  ,, , 

z  -  aj2(C(z))  +  0  aT2(C(z))(  1  -z) 


From  (5.2) 

E(oo)a-jv2(w) 
=  [  ii^ ] 

w-X2(l-(l-e+eu)aT2(w))-Xi(l-aj.i(w)) 

C(z)(z-(i-e+eu)aT2(C(z)))  ., , , 

[  — +  w  -  C(z)  ], 

z  -  aT2(^(z))  +  eaT2(C(z))(l-z) 

which,  from  (2.1)  and  (2.4) 

E(oo)a^(w)  eC(z)  (z-u)  aT2(C(z)) 

=  I  * J  [  — +  w  J 

w-(Pi+p2)wa.j.j2(w)+e(l-u)aj2(w)        z  -  ay2(C(z))  +  6aT2(C(z))(l-z) 

which,  from  (A.9), 

E(co)a^(w)  -h(z)X2e  (z-u) 

=  [  J ]  [ -^ +  w  ] 

w-(pi+p2)wa.j-j2(w)+e(l-u)a^(w)  E(oo)(l-0) 

which  simplifies  to  (5.4).* 

A-6  Proof  of  Theorem  6.1:  Using  the  argument  of  Equation  (A.  10)  one  finds  that  if  V*  is  the 
joint  distribution  of  the  duration  of  the  first  i  passes  and  N*  be  the  number  of  C2  customers  in  the 
sj'stem  at  the  end  of  the  ith  pass  if  no  other  Cj  or  Cj  customers  arrive  then 

E[u^VsV*  I  Npi=ni,  Vi  =  vj]    =  (l-0+eu)"i  aT2"i"^'(s)  e-^^l .  (A.  18) 

In  the  presence  of  future  Cj  arrivals,  joint  transform  of  the  duration  of  the  first  i  passes,  Vj*  , 
and  number  of  C2  customers  in  the  system  at  the  end  of  the  ith  pass,  Npj ,  is  found  using  the 

argument  of  (A.  15)  to  be 

E[u^P2*e-sV2*]  =  E[  u^*  e-(s  +  h  -Mobpi(s))V*  ]  (A.  19) 

where  cy3pi(s)  is  defined  in  TABLE  2.1. 
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In  the  presence  of  both  future  Cj  arrivals  and  future  C2  arrivals,  the  duration  of  the  first  i  passes, 
Vj  ,  number  of  C2  customers  in  the  system  at  the  end  of  the  ith  pass,  Npi ,  is  found  using  the 

argument  of  (A.  16),  to  be 

T^2(u,s)  =    E[u^P2  e-sV2]  =  E[u^P2*e-  (^  +  h.  -^u) V2*  ] .  (A.20) 

Combining  (A.  1 8),  (A.  19)  and  (A.20)  gives 

t52(u,s)  =  aT2(w(u,s))di(z(u,s),w(u,s)) 

where  z(u,s)  and  w(u,s)  are  defined  in  equations  (5.2)  and  (5.3). 

Using  a  similar  argument,  for  i  =  1,2,... 

i^i+i(u,s)  =  aT2(w(u,s))T3j(z(u,s),w(u,s)).  (A.21) 

If  we  uncondition  on  the  number  of  passes  then  we  find  that  the  joint  distribution  of  the  number  of 
C2  customers  in  the  system  and  the  total  elapsed  time  is 


X(u,s)  =  y  (l-e)ei-li3i(u,s)  (A.22) 

1=1 


which  equals 


=  (i-e)di(u,s)  +  e  £  (i-0)ei-idi^i(u,s) 

i=l 
which,  from  (A.21), 

CO 

=  (i-e)di(u,s)  +  eaj2(w(u,s))X  (i-e)ei-ii»i(z(u,s),w(u,s)) 

i=l 
which,  from  (A.16)  equals  (6.1)* 

A.7  Proof  of  Theorem  7.1: 

If  CJgpj(s)  is  the  Laplace  transform  of  an  busy  period  of  Cj  customers  assuming  no  C2  customer 

arrivals,  i.e., 

aBpi(r)  =  otri(r  +  X^  -  X^(5^^^{r)).  (A.23) 
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Define  (y*Bpi(s)  to  be  the  Laplace  transform  of  a  busy  period  that  begins  with  a  Cj  customers  in  the 
presence  of  both  Cj  and  C2  arrivals.  Then,  o*Bpj(s)  corresponds  to  the  duration  of  a  Tbpi 
interrupted  at  rate  X2  by  a  interruptions  having  a  duration  with  Laplace  Transform  cr*Bp2(s). 
Hence, 

«^BPi(s)  =  C^BPi^s  +  ^2  -  ^2^BP2(s))-  (A.24) 

The  busy  period  of  the  C2  customers  is  independent  of  the  order  of  service.  Therefore  it  will  be 
assumed,  as  was  done  in  Theorem  2.1,  that  the  customers  feed  back  to  the  front  of  the  line  rather 
than  the  back  of  the  line  and  therefore  have  an  effective  service  time  given  in  equation  (2.2). 

If  this  service  time  is  modified  due  to  the  interruptions  of  the  Poisson  stream  of  Cj  customers 
having  busy  periods  with  a  Laplace  transform  Ogpj(r)  then  the  modified  service  time  has  a  Laplace 

transform, 

«T2EFfW  =  "t2EFf(^  +  ^1  ■  ^i^piW)  (A.25) 

and  a  busy  period  starting  wiUi  a  C2  customer  is  given  by 

Cy*BP2(s)  =  aT2EFF(s  "*"  ^2  "  "^l^BPl^^^  ^  ■  (A.26) 

From  (A. 23)  and  (A.24),  setting  r  =  s  +  Xj  -  X2CT*Bp2(s), 

cy*Bpi(s)  =  aq<i(s  +  X2  -  ?i2a*BP2(s)  +  ^i  -  ^i<7bpi(5  "•■  ^2 '  ^2^BP2(^))) 

which  from  (A.24)  implies  (7.2). 

From  (A.25)  and  (A.26),  setting  r  =  s  +  ^.j-  ^2'^BP2(^) 

Cy*BP2(^)  =  OjjEFF^S  ■*■  ^2  "  ^2'^*BP2(^)  "^  ^1  '  ^l^BPl^^  "^  ^2  "  ^2*^BP2(^)) 

which  from  (A.26)  implies  (7.3). 

Equation  (7.1)  follows  from  the  definition  of  a3p(s),  cy*BPi(^)'  ^^^  ^bp2(^)  ^"^  ^'^•^^  ^^°"^ 
differentiations  of  (7.1),  (7.2)  and  (7.3)* 

ACKNOWLEDGMENT:  The  authors  are  grateful  to  A.  Khoo  for  tiie  programming  support 
which  was  used  to  generate  TABLE  8.1. 
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